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Abstract 

We present modeling of the incompressible viscous flows in the domain containing unconfined 
fluid and a porous medium in the case when the flow in the unconfined domain dominates. 
For such setting a rigorous derivation of the Beavers- Joseph-Saffman interface condition was 
undertaken by Jager and Mikelic [S1AM J. Appl. Math. 60 (2000), p.1111-1127] using the 
homogenization method. So far the interface law for the pressure was conceived and confirmed 
only numerically. In this article we derive the Beavers and Joseph law for a general body force 
by estimating the pressure field approximation. Different than in the Poiseuille flow case, the 
velocity approximation is not divergence-free and the precise pressure estimation is essential. 
This new estimate allows us to justify rigorously the pressure jump condition using the Navier 
boundary layer, already used to calculate the constant in the law by Beavers and Joseph. Finally, 
our results confirm that the position of the interface influences the solution only at the order of 
physical permeability and therefore the choice of this position does not pose problems. 
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1 Introduction 



Slow viscous and incompressible simultaneous flow through an unconfined region and a porous 
medium occurs in a wide range of industrial processes and natural phenomena. One of the classical 
problems is finding effective boundary conditions at a naturally permeable wall, i.e., at the surface 
which separates a channel flow and a porous medium. 

The effective laminar incompressible and viscous flow through a porous medium can be described 
using the Darcy's law. The unconfined fluid flow in the channel is governed by the Stokes system, 
or by the Navier-Stokes system if the inertia effects in the free fluid are important. To model 
the coupling of both processes, it is necessary to put together two systems of partial differential 
equations: the second order system for the velocity and the first order equation for the pressure, 

-/xAu + Vp = / (1) 
div u = 0, (2) 

in the unconfined fluid region, and the scalar second order equation for the pressure and the first 
order system for the seepage velocity, 

-»v F =K(f-Vp F ) (3) 
div v F = 0, (4) 

in the porous medium. 

The orders of the corresponding differential operators are different and it is not clear what 
conditions it is necessary to impose at the interface between the free fluid and the porous part of 
the domain. One coupling condition is based on the continuity of the normal mass flux. However, it 
is not enough for determination of the effective flow and it is necessary to specify more conditions. 

Several laws of fluid dynamics in porous media were derived using homogenization. The most 
notable example is the Darcy's law, being the effective equation for one phase flow through a rigid 
porous medium. Its formal derivation using the 2-scale expansion goes back to the classical paper 
by Ene and Sanchez-Palencia [8]. This derivation was made mathematically rigorous by Tartar in 
reference [24] . For the detailed proof in the case of a periodic porous medium we refer to the review 
papers by Allaire [1], and by Mikelic [18] and for a random statistically homogeneous porous medium 
to the paper of Beliaev and Kozlov [3] . 

As in the derivation of Darcy's law, we would like to apply the homogenization technique to 
find the effective interface laws. However, the assumption of statistical homogeneity of the domain, 
which is necessary for the homogenization approach, is not valid close to the interface. Consequently, 
deviations from the Darcy's law are expected in the thin layers near the interfaces. Furthermore, 
presence of such interfaces can significantly change the structure of the model coefficients and lead 
to different effective constitutive laws for the flow. 

It was experimentally found by Beavers and Joseph in [2] that the jump of the tangential compo- 
nent of the effective velocity at the interface is proportional to the shear stress originating from the 
free fluid. This law was justified at a physics level of rigor by Saffman in [21], where it was observed 
that the seepage velocity contribution could be neglected leading to the law in the form 

Vk^ =av T + 0(k), (5) 

where a is a dimensionless parameter depending on the geometrical structure of the porous medium, 
e is the characteristic pore size, and k is the scalar permeability, v denotes the unit normal vector 
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at the interface and v T is the slip velocity of the free fluid in the channel. Saffman's modification of 
the law by Beavers and Joseph has been widely accepted. 

As an alternative to (5), the continuity of the effective pressure was suggested by Ene and 
Sanchez-Palencia in [8]. While this interface law is acceptable from modeling point of view, it 
should be noted that the well-posedness of the averaged problem is not clear. 

The law (5) was rigorously justified by Jager and Mikelic in [12]. Numerical calculations of the 
boundary layers for the experimental conditions of Beavers and Joseph are presented in [13]. They 
indicate appearance of a pressure jump at the interface. These issues were heuristically discussed in 
[14]- 

In the experiment by Beavers and Joseph only the flows tangential to a naturally permeable 
wall (a porous bed) were considered. In general, the situation is much more complicated and 
many types of interfacial conditions have been proposed, such as continuous tangential velocity with 
discontinuous tangential shear stress introduced in [20] by Ochoa-Tapia and Whitakcr, or continuous 
tangential velocity and tangential shear stress in reference [19] by Neale and Nader, or discontinuous 
tangential velocity and tangential shear stress from [5] by Cieszko and Kubik. In particular, in 
[20] the continuity of the velocity and the continuity of the " modified" normal stress were obtained 
at the interface using volume averaging. In order to perform the averaging it was necessary to 
assume the Brinkman's flow in the porous part and a transition layer between the two domains. 
Numerical study of the hydrodynamic boundary condition at the interface between a porous and 
plain medium was performed by Sahraoui and Kaviany [22]. Numerical implementation of the 
effective interface couplings was presented in [7] and in [10]. Nevertheless, determination of the 
practical and relevant first-order interface conditions between the pure fluid and the porous matrix 
remains an open question that could be treated using the technique developed in reference [11]. 

This paper is a continuation of works [12] and [13] and constitutes a step forward in the de- 
velopment of the rigorous approach to model effective interface laws for the transport phenomena 
between an unconfincd fluid and a porous medium. We depart beyond justification of the law (5) 
developed in [12] and undertake a rigorous derivation of the interface laws for the viscous flow in 
a long channel in contact with a porous bed. The macroscopic model derived links pressure jump 
with the shear stress of the unconfined fluid at the interface, an effect which was predicted based 
on numerical simulations in reference [13]. Derivation of the law of Beavers and Joseph is based 
on the procedures proposed in [12] and discussed in [14]; however it is nontrivially adjusted to the 
new setting involving a general body force. We consider a general situation when the flow in the 
unconfined region dominates. Nevertheless, even if the flow is much less important in the porous 
part, the pressures are of the same order of magnitude. Hence finding and justifying the interface 
law for the pressure is of fundamental interest. 

The review paper [14] was concluded with the sentence "Proving the error estimate for the 
pressure approximation in the porous bed f2| remains an open problem" . We solve this problem 
and present a mathematically rigorous derivation of the pressure jump interface law, which is the 
next order correction of the Beavers- Joseph law. We obtain the effective equations heuristically and 
then rigorously justify them. Combination of homogenization and boundary layer approaches is 
used to achieve this end. Study of such complex flows leads to artificial compressibility effects in the 
upscaling process. In this paper we develop the required estimate of the pressure. Our main results 
are the following: 

1. Confirmation of Saffman's form of the law by Beavers and Joseph in the more general setting 

u{ ff = -ecf-^+0{e% (6) 
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where u e ^ is the average over the characteristic pore opening at the naturally permeable wall. 
Physical permeability is given by k — k E = e 2 K and the constant in (6) is proportional to 
The error is of order k e , as remarked by Saffman in [21]. It is important to point out that 
the parameter a from expression (5) is determined taking into account the auxiliary problems, 

which we formulate later in (105)-(108) and (111), and that it is given by a — — > 0. 

2. Interface between the unconfincd flow and the porous bed is an artificial mathematical bound- 
ary and it can be chosen in a layer having the pore size thickness. We show that a perturbation 
of the interface position of the order 0{e) implies a perturbation in the solution of 0(e 2 ). Con- 
sequently, it influences the result only at the next order of the asymptotic expansion. 

3. We obtain a uniform bound on the pressure approximation. Furthermore, we prove that there 
is a jump of the effective pressure on the interface and that it is proportional to the free fluid 
flow shear at the boundary. The proportionality constant is calculated from the boundary 
layer problem (105)-(108). Homogenization leads to the discontinuity of the effective pressure 
field at the interface, which differs from the pressure interface continuity law proposed in 
reference [8] . If the boundary layer pressure is neglected, the pressure in the neighborhood of 
the interface is poorly approximated. 

Here, we remark that some classes of problems, like infiltration into the porous medium, are charac- 
terized by the velocity field of the same order in both domains. Such situation requests much larger 
body force in the porous part than in unconfincd. Some situations of this kind were considered in 
[11]. In this paper, the body force is of order 0(1) in both domains. 

The paper is organized as follows. In Section 2 we formulate the problem and main results. 
Section 3 is devoted to the proof of the results. We conclude the paper with two short appendices 
recalling the notion of very weak solutions and definition and properties of the Navier boundary 
layer. 

2 Statement of the problem and of the results 
2.1 Definition of the geometry 

Let L, h and H be positive real numbers. We consider a two dimensional periodic porous medium 
Sl 2 = (0, L) x (— H, 0) with a periodic arrangement of the pores. The formal description goes along 
the following lines: 

First, we define the geometrical structure inside the unit cell Y = (0, l) 2 . Let Y 8 (the solid part) 
be a closed strictly included subset of Y, and Yp = Y\Y S (the fluid part). Now we make a periodic 
repetition of Y s all over R 2 and set Y^ =Y s + k, k eZ 2 . Obviously, the resulting set E s = (Jfcgz 2 
is a closed subset of R 2 and Ep = R 2 \i? s in an open set in R 2 . We suppose that Y s has a boundary 
of class C 0,1 , which is locally located on one side of their boundary. Obviously, Ep is connected and 
E s is not. 

Now we notice that ^2 is covered with a regular mesh of size e, each cell being a cube Yf, with 
1 < i < N(e) = |il2|e~ 2 [l + o(l)]. Each cube Yf is homeomorphic to Y, by linear homeomorphism 
Ilf , being composed of translation and a homothety of ratio 1/e. 

We define F|. = (ITf)" 1 ^) and Yp. = (nf)" 1 (y F ). For sufficiently small e > we 
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Figure 1: The geometry 



consider the set T e = {k e Z 2 \Y§ k C f^} and define 

°* = U y s fc , s E = dO £ , fi§ = 2 \0 £ = Q 2 n £ £ F 

fe£T e 

Obviously, <9£1| = <9£l2 U S e . The domains O e and f2| represent, respectively, the solid and fluid 
parts of the porous medium Q. For simplicity, we suppose L/e, H/e, h/e € N. 

We set E = (0,L) x {0}, fti = (0,L) x (0,/i) and SI = (0,i) x (-H,h). Furthermore, let 
ST =fi|UEUfii. 

A very important property of the porous media is the following variant of Poincare's inequality: 
Lemma 1. (see e.g. [23]) Let ip e V(fi|) = G i? 1 ^!) I<P = on Se )- Then > u holds 

IMU=»(e) < C , £ 1 / 2 ||V x ^|| L 2 (n|) 2, (7) 
IMU 2 (fi|) <Ce\\V x ip\\ L 2 (ni) 2. (8) 

2.2 The microscopic equations 

Having defined the geometrical structure of the porous medium, we precise the flow problem. Here 
we consider the slow viscous incompressible flow of a single fluid through a porous medium. We 
suppose the no-slip condition at the boundaries of the pores (i.e., a rigid porous medium). Then, 
we describe it by the following non-dimensional steady Stokes system in Sl e (the fluid part of the 
porous medium fi): 



-Av e + Vp £ = f 
divv £ = in ST, 

v £ = on <9ST \ \{xi = 0} U {xi = L} 



in sr 

p e dx = 0, 



/ 

J Six 



(9) 
(10) 



{v e ,p e } is L — periodic in x\. (11) 



Here the non-dimensional f stands for the effects of external forces or an injection at the boundary 
or a given pressure drop, and it corresponds to the physical forcing term multiplied by the ratio 
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between Reynolds' number and Froude's number squared. v £ denotes the non-dimensional velocity 
and;/ is the non-dimensional pressure. The non-constant force / corresponds, e.g., to a non-constant 
pressure drop or to injection profiles which are not parabolic. 
Let 

W e = {z e H 1 ^) 2 , z = on dil e \ ^{x 1 — 0}U{x 1 = L}^j and z is L -periodic in xj. (12) 

The variational form of the problem (9)-(ll) reads: 

Find v £ e W e , div v £ = in £ and p £ <= L 2 (D, e ) such that 

/ Vv £ V^dx- / p e div (pdx= I fipdx V(peW e . (13) 
Jn e Jn E Jn e 

Then for f e C°°(f2) 2 , the elementary elliptic variational theory gives the existence of the unique 
velocity field v e <G W e , div v e = in f2 e , which solves (13) for every <p e W £ , div ip = in ft 6 . The 
construction of the pressure field goes through De Rham's theorem (see e.g. book [25]). 



2.3 Main result 

We start by introducing the effective problems in fii (the unconfined fluid part) and Q.2- 
Find a velocity field u° and a pressure field p e ^ such that 

-Au e// + Vp eff = f in Sl u (14) 



divu e// = inOi, / p eff dx = 7 



(15) 



u eif = on (0,L) x {h}; u eff and p eif are L - periodic in x\, (16) 
u e 2 ff = Q and u\ ff +eC\ 1 ^- = on £. (17) 

We note that the second boundary condition in (17) is the law by Beavers and Joseph from [2]. 
The constant C± l is strictly negative and calculated through (111), from the viscous boundary layer 
described in Appendix 2. 

Problem (14)-(17) has a unique solution, which in the case of Poiseuille flows (i.e. when f = 

Pb-Po n , 
e ) reads 



eft _ [ Pb-Po ( _ eC\ l h 
pms \ 2L { 2 h-eC\ 



c i, u (x2-h),QJ for < x 2 < h; p eff = for < x 1 < L. (18) 
The effective mass flow rate through the channel is then 

M eff = [ u\ n dx, (19) 

which for the Poiseuille flow reads M e JJ s = - Pb ~ P " h 3 k ~ . (20) 

LZ ft £L/^ 
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Theorem 2. Let us suppose f G C°°(VL) 2 and L-periodic with respect to X\. For {\r £ ,p £ } given by 
(9)-(ll) and{u e ff,p e ff} by (U)-(17). It holds 

\\v £ - u e ^|| i2(ni)2 + \M £ - M e "\ < Ce 3 ' 2 (21) 

\\v £ - u e "|| ffl , w + Hp 6 -p e ^|| L i (ni) + ||V(v e - u e ")IUw + 

IIN^V^ - u e ")IUw + HN V V -P eff )\\L H n^ < Ce, (22) 

with M e ff defined in (19). 

Next, we study the situation in the porous medium Q, 2 . 

Theorem 3. Let the permeability tensor K be given by (83). The effective porous media pressure 
p° is the L— periodic in X\ function satisfying 

div (^K(f(x) - Wp°)j =0 in fl 2 (23) 

f=p eff + CS^-{ Xl ,0) on S; K(t (x) - Vp°)\ {x2= _ H} ■ e 2 = 0, (24) 

with u e ff being the solution to the problem (14)- (17) and being the pressure stabilization constant 
defined by (113). In addition we have 



1 

e 

p £ -p°^0 strongly in L 2 (fl 2 ), as e -> 0; (26) 



2 v £ -K(f - Vp°) weakly in L 2 ((0,L) x (-H,-5)) 2 , as e -S- 0, > 0; (25) 



M-p efI \\ H -u^)<C^. (27) 



Remark 4. // we include the vicinity of £ the velocity v e has to be corrected by a boundary layer 
term f3 U ' £ (x) — e(3 bl (x/e), defined through (105)-(108), and the convergence result (25) reads 

_( v e + / gw, e __l_( a . lj0 )) -K(f(x)-\7p°)^0 weakly in L 2 (n 2 ) 2 , as e -> 0. (28) 

Remark 5. Let fl ae = (0, L) x (ae, ft) for a < and let { u a > e f f , p a > e f f } be a solution for (14)-(17) 
in £l a e, with (17) replaced by 

f)„. a > e ff 

u a ^ n = and u a { ett + eCl' bl ^ - on S a = (0, b) x ae. (29) 

Problem (14)-(16), (29) has a unique smooth solution {u a ' e ^ ,p a ' e ^} , its derivatives are bounded 
independently of e and, by (117), C°^ hl = C\ l — a. Then a simple calculation gives 

= u a ^ S (xi , eo) + eC?' bl 9v ^- (an , ea) = u a { eff (an , 0) + eC? (an , 0) + 

f F \2 a2 a,e// o2 a,eff 

^-(^5-( a; i>^) + £L ^-(^6)). /° r a,6e(0, £ a). 

Therefore, a perturbation of the interface position for an 0(e) implies a perturbation in the solution 
of 0(e 2 ) in H k (Qi). Consequently, there is a freedom in fixing position o/S. It influences the result 
only at the next order of the asymptotic expansion. 
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The physical permeability K p h ys is proportional to e 2 . Our result on the influence of the interface 
position on the effective slip is in agreement with the observation of Kaviany in [15], pages 79-83. 
In fact, it has been noticed by Larson and Higdon in [16] that changes of 0(1) in the slip coefficients 
are possible, after the change of order 0(y/K p i iys ) of the interface position. Therefore, the exact 
position of 12 does not pose problems, since it influences the solution only at order 0(K p h ys ). 



3 Law by Beavers and Joseph 

In this section we extend the justification of the law (5) from [12] to the case with a general body 
force. Our boundary conditions are simpler from those of the experiment from [2] and we consider the 
2D Stokes system. The Beavers and Joseph setting could be reduced to our setting if £1 is sufficiently 
long in x\ direction. Then we may assume the periodic boundary conditions at inlet/outlet boundary 

and the flow is governed by a force coming from the pressure drop and is equal to ^——^—e 1 . We 

o 

assume a non-constant force, which can describe a larger class of the problems. 
3.1 The impermeable interface approximation 

Intuitively, the main flow is in the unconfined domain f2i . Following the approach from [12] we 
study the problem 

-Av° + Vp°=f infii, (30) 
div v° = in n u (31) 

v° = on dfij \ ^{x-l = 0} U {x-l = L}^j , (32) 

{v°,p } is L — periodic in x\ (33) 

Problem (30)-(33) has a unique solution {v°,p } e H 1 ^) 2 x Lg(fii) (see e.g. book [25]). In fact 
this solution is C°° for f e C°°. Therefore, for the lowest order approximation {v°,p } we impose 
on the interface the no-slip condition 

v° = on S. (34) 
We observe that in the Beavers and Joseph setting f = — Pb Po e 1 and the unique solution for 

Li 

this problem in _ff 1 (r2 1 ) 2 x Lq(SIi) is the classic Poiseuille flow in fli, satisfying the no-slip conditions 
at 12. It is given by 

v° = ^ Pfc 2 ~ P0 x 2 (x 2 - ft), 0^ for < x 2 < h; p° = for < x x < L (35) 

(see [12] and [14] for further details). 

We extend v° to r2 2 by setting v° = for —H < x 2 < 0. For p° we use a smooth extension to ^2, 
p , which we shall precise. The question is in which sense this solution approximates the solution 
{v e ,p £ } of the original problem (9)-(ll). 
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Direct consequence of the weak formulation (13) is that the difference v £ — v° satisfies the 
following variational equation: 

f V(v £ -v°)V<^ dx- [ (p £ -p°) div<p= [ l^^j dS- [ [p a }f 2 dS+ [ (f-Vp°)<p dx, V<p e W £ . 

(36) 

Taking <p = v £ — v° in (36) and applying Lemma 1 leads to the following result, proved in [12]: 

Proposition 6. Let {v £ ,p £ } be the solution for (9)-(ll) and {v°,p } defined by (30)-(33). Then, 
it holds 

^||V(v £ - v°)\\ L 2 {nc)i + -L||v £ || L 2 (a|)2 + ||v e || L 2 (E) < Ce (37) 

Furthermore, using estimate (37) and the notion of very weak solutions for the Stokes system in 
ill, introduced in [6] (see also Appendix 1), we conclude the following additional estimates: 

Corollary 7. (see [12]) Let {v £ ,p £ } be the solution for (9)-(ll) and {v°,p } defined by (30)- (33). 
Then, it holds 

yj~e\\p £ - p°\\ L 2 {ni) + ||v £ - «°|| La(ni)a < Ce. (38) 

This provides the uniform a priori estimates for {v £ ,p £ }. Moreover, we have found that the 
viscous flow in f2i corresponding to an impermeable wall is an 0(e) L 2 -approximation for v e . Beavers 
and Joseph's law should correspond to the next order velocity correction. Since the Darcy velocity 
is of order 0(e 2 ) we justify Saffman's version of the law. 

3.2 Justification of the law by Beavers and Joseph 

At the interface S the approximation from Subsection 3.1 leads to the shear stress jump equal to 

dv° 

— — Mj> Contrary to the pressure difference, which could be easily set to zero by the appropriate 
0x2 

choice of p°, the shear stress jump requires construction of the corresponding boundary layer. For 

the intuitive argument how to obtain the shear stress jump correction using the natural stretching 

x 

variable y = — , we refer to the paper [14], page 503. In the present paper we present the rigorous 

e 

construction, based on the Navier boundary layer and following the scheme originally used in [12]. 
Let {P bl ,u) bl } be the boundary layer given by (105)-(108). 
Now we set 

P bl ' s {x)=sp bl (^) and J l ' £ (x) = w fc '(|), xen £ , (39) 
f3 ' e is extended by zero to SI \ il £ . Let H be Heaviside's function. Then for every q > 1 we obtain 

^\\f3 u ' £ - e(C b 1 l ,0)H(x 2 )\\ Lq(nr + \\u bl ' £ - C bl H(x 2 )\\ L « m + ||V/? w > £ |U, (ni uEufi 2 )4 - Ce 1 ^. (40) 

Hence, our boundary layer is not concentrated around the interface and there are some stabilization 
constants. We will see that these constants are closely linked to our effective interface law. 
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As in [11] stabilization of f3 bl ' £ towards a nonzero constant velocity e(C bl , 0), at the upper bound- 
ary, generates a countcrflow. It is given by the following Stokes system in Q,\: 

-Az CT + Vp°" = infii, (41) 

divz CT =0 infii, (42) 
dv° 

z ff = on {x 2 = h} andz' 7 = — J-he 1 on {x 2 = 0}, (43) 

ox 2 

{z° ,p a } is L — periodic in x\. (44) 

In the setting of the experiment by Beavers and Joseph, z CT was proportional to the two dimensional 

Couettc flow d = (1 - ^)e 1 . 

h 

Now, after [11], we expected that the approximation for the velocity reads 

f)v° 

v e = v° - - £ (Cf ,0))g| E - eCfz* + 0(e 2 ), (45) 

Concerning the pressure, there are additional complications due to the stabilization of the bound- 
ed 

ary layer pressure to Cjj , when y 2 — > +oo. Consequently, uj bl ' 6 — H (x 2 )C^ 0~~He is small in fli and 

we should take into account the pressure stabilization effect. 

At the flat interface S, the normal component of the normal stress reduces to the pressure field. 
Subtraction of the stabilization pressure constant at infinity leads to the pressure jump on S: 

[p E h^P°(x 1 ,+0)-p°(x 1 ,-0) = -C b J^\j: + O(e) for Xl G(0,L). (46) 

Therefore, the pressure approximation is 

p%x) = p°H(x 2 ) +p°H(-x 2 ) - (w bl > s (x) - H(x 2 )C b J)^k ~ sC bl p°H(x 2 ) + 0(e). (47) 

Following the ideas from [11], these heuristic calculations could be made rigorous. Let us define 
the errors in velocity and in the pressure: 

W{x) = v e -v + (f3 bl ' e -eC bl e 1 H{x 2 ))^\i;+eC b 1 l 2,' 7 (48) 

V s (x) =p s -p°H(x 2 ) -p°H(-x 2 ) + (cu bL %x) - H(x 2 )C b J)^-\x + sC bl p° H (x 2 ) . (49) 

Remark 8. Rigorous argument, showing that U e is of order 0(e 2 ), allows justifying Saffman's 
modification of the Beavers and Joseph law (5): On the interface £ we obtain 

dv\. dv9. d(5 bl . . , v\ nhl , , nS dv9. ^. . 



dx 2 IZj dx 2 ^ dy 2 x ^ 6 " s ^ l dx 



2 



After averaging over £ with respect to yi, we obtain the Saffman version of the law by Beavers and 
Joseph 

u e / f = -eC bld $— on E, (50) 
ox 2 

where is the average of vf over the characteristic pore opening at the naturally permeable wall. 
The higher order terms are neglected. 
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For simplicity we denote 



/ \ 9V U 

■bl,s r-nbl^l EJY~ w 



Then, the variational equation for (8 ' £ — eC\ e H(x 2 )) — — Is reads 

ox 2 

f v((B u - £ - eC'l l e 1 H(x 2 ))<j° 12 \ :X7<pdx- f a? 2 (w W ' £ (x) - H(x 2 )C b J) div <p dx = 
- [ Vl a° 12 dS- [ Ct l V2a 12 dS- [ V ( Aa° 12 ((3 bl ' e - eC bl 5 u H{x 2 )) Vi - 

d Xi a° 12 (w bl ' e - eC2)<Pi - 2(8 bl ' £ - eC\ l 5 u H(x 2 )) div (^W° 2 )) dx, e W e . (51) 
Next, the variational form of (41)-(44) reads 



Vz, a :V(p dx- ( p a div tp dx = -[ (-ip 2 p a + (p-^) dS, Vtp eW £ . (52) 



Now we are ready to write the variational equation for {U e ,V e } and obtain the higher order 
error estimates as in [12]. Nevertheless, contrary to [12], U e is not divergence free anymore and we 
need more effort to control V e . 

Theorem 9. LetU e be defined by (48) andV e by (49). Let p° be a smooth function satisfying the 
interface condition (46). Then, the following estimates hold 

e||V7> e ||H-M^) +4^U £ \\mn-Y + \\U e \\mn^ + £ 1/2 ||^|| L 2 (E)2 < Ce 2 (53) 

Proof. First we remark that for y 2 > the mean with respect to y\ of uj u (y) — C% is zero. Conse- 
quently, the problem 

dir bl f 1 

= u bl (y) - C b J Vyi G (0, 1); n b J is 1-pcriodic ; / n b J( yi ,y 2 ) d Vl = 0, (54) 

oyi Jo 

has a unique smooth solution. 

Next by subtracting (51) and (52) from (36) we obtain 

/ VIA 6 : Vtp dx - / V £ div tp dx = 

£ I {-V2P° + <fi ■ dS + [ (f - Vp°V dx- [ V ( Aa° 12 (/3- L£ - sC bl 6 li H(x 2 ))<p i dx- 
Jt, ox 2 J s J U e i -r J \ 

f d Xl ol 2 J*> e tpidx- [ sw b J( X )(i Pl d 2 Xi a 12 + d Xl <p 1 d Xl a° 12 )dx+ 

2 f (Bf e - eC bl H(x 2 ))( Vl d 2 Xi a° 12 + d Xl ^d Xl a° 12 ) dx, e W £ , (55) 



&wU e = {I3 b j> £ - sCl l H{x 2 ))-^o-° l2 in fl E . (56) 
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From (55) we find out that 

| / VW:V<pdx- [ V e div <pdx\ <Ce 3 / 2 ||V^|| L 2 (n£) 4+C £ ||f-Vp |U2 (n2)2 ||V^|U2 (ni) 4 

(57) 

and || div W £ || L 2 (OE) 2 < Ce 3/2 . (58) 

The size of div U e does not allow us to obtain the appropriate estimate and we should diminish it 
further. 

Let Q w be given by (118)-(120). Furthermore let Q bl ' E (x) = e 2 Q hl (x/e) and let w Q be defined 

by 

f Aw Q - Vp Q = in 



div w<3 = j^-r [ 4-^12 dS = in O i; 



wQ = --j—v^e 2 on E, w Q = on {x 2 = h}; 



(59) 



, {w^jp 5 } is L-pcriodic in x\. 
We introduce the following error functions, where the compressibility effects are reduced to the next 



order: 



Wix) = W {x) + Q w ' e (aO-^-a? 2 + e 2 H(x 2 )( f {C\ l H{y 2 ) - ^(y)) dy)w« , (60) 

" x i Jz BL 

T%x)=V^x,t) + s 2 H(x 2 )([ {Cl l H{y 2 )-${y))dy)pQ, (61) 

(62) 



div^ £ = -Qf' £ (x) -^<7° 2 in ft £ . 



Then Uq e W E and || div Wo||l 2 (o») 4 < Ce 5/2 . Next, we construct a function $ 1,e e H 1 ^) 2 sucn 
that 

div &<* = -Q\ l ' £ (x) ^*° 12 in Q u ; 



$ 1,£ = 1^7 / Qi' e (x)-^a% 2 dx on E, $ l£ = on {x 2 = h}, 
1^1 Jn 1 dx 1 

$ 1 ' e is L-periodic in X\. 



(63) 



We note that ||$ 1,£ ||ffi(Oi) 2 < C^ 2 - Next wc extend Q H ' £ by zero to the rigid part of the porous 
medium and choose a function $ 2 ' £ <G i? 1 (ri 2 ) 2 such that 



div $ 2 < £ = -Qf e {x) -^a° 12 in 2 , 



<J>2,£ 



dx 2 

/ Qr e (x)—^a^ dx on E, $ 2 < £ = on {x 2 = -H}, 
JQ.2 dx 1 



(64) 



<f> 2,£ is L-periodic in x\. 
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We note that $ ll£ = $ 2e on £ and ||$ 2 ' e || H i (ni) 2 < Ce 2 . Let I 2 = {z£ H 1 ^) 2 , z = on {zi = 
L} and z is L — periodic in xi} and Xf = {z G X 2 , z = on 9Q 2 \ df2 2 }- In the seminal 
paper [24] Tartar constructed a continuous linear restriction operator operator R e G C(X 2 , X 2 ), 
such that 

div (R E <p) = div ip + V] ^ x ye / div 95 cfe, V99 G X 2 (65) 

H^vlU^ns)" < C{£||V(^|| L 2 (n2) 4 + M L 2 in2)2 }, \/<p g X 2 (66) 
C 

||V(i? e (^)|| L 2 ( ^ ) 4 < j{e||Vv||£2(n a )« + IMU'(n a )'}. G X 2 . (67) 

Furthermore, ip = R E ip on E. For more details we refer also to [1] and [18]. This construction allows 
us to work with the divergence free velocity error function U given by 

U £ = Ul - H(x 2 )^' £ - H{-x 2 )R e ^ £ (68) 

Now we write the analogue of the variational equation (55) for {if , P e } and, since | |Vi? e $ 2 ' e | 
< Ce. We find out that the leading order force term is of the same order as in the estimate (57). 
Now we test the analogue of variational equation (55) for {U ,V £ } with ip = U to obtain 

||W £ || L 2 (0e) 4 < Ce. (69) 

We remark that if differs from W for 0(e 2 ) in L 2 -norm and for O(e) in iJ 1 -norm . Therefore (69) 
gives us the middle part of the estimate (53). In what concerns the L 2 {Ti) norm of U e , it follows by 
using (7). Remaining pressure estimate follows easily from the weak formulation and the estimates 
onW £ . ' □ 

Next we use Theorem 9 and the results on the Stokes system with L 2 — boundary values from 
[9] and [4] to conclude the following result: 

Corollary 10. Let U e be defined by (48) and V s by (49). Let p a be a smooth function satisfying 
the interface condition (46). Then, the following estimate holds 

V~e\\V £ \\ LHni) + ||« e ||*i/ w < Ce 3 / 2 . (70) 

Now we introduce the effective flow equations in £li through the boundary value problem (14)- 
(17), containing the slip condition of Beavers and Joseph. Since our expansion is performed using the 
solution {v°,p } of the problem (30)-(33), we need to know the relationship between the solutions 
to these two boundary value problems. 

Proposition 11. Let f G C°°(Sli) 2 and L-periodic in x\. Let {u e ff ,p e ^} be the solution of the 
problem (U)-(17), {v°,p } of the problem (30)- (33) and {z a ,p a } of the problem (41)-(44). Then 
we have 



u 



eff 



/ ) ||ff*(n 1 )» + lb e// -P°llff*-i(n 1 ) <Cs, VfcGN; (71) 
l|u e// - v° + eC?z°\\ HHsil) 2 + \\p eff -p° + eC b M\\ Hk - 1(ni) < Ce 2 , Vfc G N. (72) 

Proof. The elliptic regularity for the Stokes operator (see e.g. [25]) gives C°° regularity for the 
functions {u e ^ ,p e ^}, {v°,p } and {z a ,p a }. It is easy to see that {u e ^ ,p e ^} is bounded in 
H k {Qi) i , independently of e, for every integer k. 
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Let U = u e t f — v° and P = p ef 1 — p°. Then for every ip e V = {<p g H 1 (£l 1 ) 2 \ tp is L-periodic 
in xi, <p = on {.X2 = h}, Lp 2 = on £} we obtain 



/ VU :\7tp dx- [ P div p dx / Uipi dS = - [ u\^p x dS. 

Using p = U as a test function yields 



(73) 



|U||j/i( ni )2 + -^=||C/i|| L 2 (E) < CVe, ^ 



1^11^(00 < Cy/l. 
Differentiating the equations with respect to x\ leads to the estimate 

,dU.. 1 ..dU x .. _ / _ 



(75) 



9[/ 

Since — — = on E, we have for the velocity trace U g i7 1 (S) 2 and its norm is smaller than Ce. 
Using [9] and [4] we obtain that 

l|U||ij3/2 (ai )2 + ||-P||jyi/2( ni ) < Ce. (76) 

After bootstrapping, we conclude that the estimate (71) holds true. 

Using corrections U 1 = u eff - v° + eCf z CT and P 1 = p eff - p° + eCfp" , for every tp e V = 
{<p g i? 1 (r2i) 2 | ^3 is L-periodic in X\, ip = on {x 2 = /i}, <p>2 = on £} we obtain 

/ VU 1 : Vif dx — f P 1 div p dx - / C/^i dS = e [ gtpi dS, (77) 

9z CT — 

where g = — C 1 ^~^|s € C°°(S) is uniformly bounded with respect to s. Repeating the argument 
used in the first part of the proof to {U 1 , P 1 } yields the estimate (72). □ 
Proof, (of Theorem 2 ) We remark that on £ 

f)v° 

v £ - u eff =W - {/3 bL£ - e (cf ,0))|^(a:i,0). (78) 

Now Theorem 9, Corollary 10 and Propositions 15 and 16 from the Appendix 1 imply the desired 
result. □ 



3.3 Justification of the interface pressure jump law and the effective equa- 
tions in the porous medium 

We have already seen that, after extension by zero to the rigid part, the velocity U e satisfies the 
a priori estimates (53), (70), with f2 e replaced by O. Furthermore, it would be more comfortable 
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to work with the pressure field V e defined on fl. Following the approach from [17], we define the 
pressure extension V £ by 

- e J V e in fl £ 

^ = \ Wl frt ve in the Y ^ for each ( ~ 79 ' ) 

where Yp. is the fluid part of the cell Yf. Note that the solid part of the porous medium is a union 
of all Yg, . Then, following Tartar's results from [24] we have 

< VP £ , ip >o=< W e ,R eV >qs, \/<p e ff 1 ^) 2 , 

where 

~ , . / <p(a:), for xerJiUS; . . 

= 1 RMx), for x e fig. ( 8 °) 

Using the estimate (53) and properties (65)- (67) of the restriction operator f? e , we arrive at 

Corollary 12. fa priori estimate for the pressure field in 0.2)- Let V s be defined by (79). Then it 
satisfies the estimates 

\\VP e \\ w ><C and \\f £ \\ LH n 2) <C, (81) 
where W = {z e H 1 (il 2 ) 2 : z = on {2:2 = — if} U {x 2 = 0}, z is L — periodic}. 

We remark that in ^2 we have strong L 2 -compactness of the family {V s }. From the properties 
of Tartar's restriction operator (see [24] or [1]) it follows: 

Lemma 13. The sequence {V s } is strongly relatively compact in L 2 {n 2 ). 

Following the homogenization derivation of the Darcy law from [8], [24], [1] or [18], we consider 
the following auxiliary problems in Yp: 

For 1 < i < 2, find {w J , tt 1 } e H^ er {Y F ) 2 x L 2 (Y F ), J Yp ir l (y) dy = 0, such that 

-A y w i (y) + V y n i (y) = e l in Y F 

div y w l (y) = in Y F (82) 

w i {y)=0 on (dY F \dY) 

Obviously, these problems always admit unique solutions. Let us introduce the permeability 
matrix K by 

Kij = [ V a w l : V y w j dy= [ w) d Vl 1 < i,j < 2. (83) 

Then after [23], permeability tensor K is symmetric and positive definite. Consequently, the drag 
tensor K^ 1 is also positive definite. 

Proof. (Proof of Theorem 3) Let the function p° be the solution for the boundary value problem 

div (^K(f(x) - Wp°)^j = in n 2 (84) 
p°=p + C b Jo-° 12 (xi) on E; K(f (x) - \7p°)\ {x2= _ H} ■ e 2 = 0. (85) 
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We take as test function in (55) ip{x)ij){y) , with ip e Qj^f^) and ip € Hp er (Yp) 2 , div y ip = 0. 
Then after passing to the subsequence 

^ -> W imp (ar, y), -»■ V^ mp (:c, y) and V e -> 7>"™ p (x) 

and we have 

/ / V y U tmp :\7yipipdydx- ( [ V tmp (x)ip(y)\7 x ^(x) dydx = [ [ (f-Vp°)il>(y)(p(x) dydx, 
Jn 2 Jy f Jn 2 Jy f Jn 2 Jy f 

(86) 

implying 

2 d(v° + V mp ) 

w mp (x, y ) = J2 wJ (y)(M x ) — -) ( a - c -) in ^- ( 87 ) 

j=l 3 

Consequently, we obtain p° + V" np G H 1 ^). 

By Corollary 10 it holds that e _1 VW e strongly in L 2 ^) 2 . Next taking tp e C£°(n) and 
using a priori estimates (53), the variational equation (55) yields a generalized form of (86) leading 
to 

V imp = on S. (88) 

Averaging div U e in Q2 results in 

div (K{i{x)- \7(V mp +p°)^j =0 in fi 2 . 
Hence the function -p jm P +p° is L— periodic in x\ and satisfies 

div (k{{{x) - \7(P mp + f))j =0 in fi 2) (89) 

P imp +p 4, =p°+^a? s ( I1 )onS; ^(f(x) - V{V imp + P°))\ {x2= - H} ■ e 2 = 0, (90) 

and we have r p imp — 0. Let p° be the solution to the problem (23)-(24). Using Proposition 11 we 
find out that p° and p° differ for Ce in any H k (fl 2 ), k e N. Hence we have established (25)-(26). It 
remains to prove the last stated result i.e. 

W-p°\\H-V*<n<Cy/i. (91) 

We use the variational equation (55) with test function having support in fii and Corollary 10 
to obtain 

|| div (VWI-^e 2 -2^ £ ^e 1 )|| i2(0l) + ||VW!-P £ e 2 -2^ e ^e 1 || L2(fil)2 <Ce. (92) 
Estimate (92) implies the following estimate for the trace 



dU. 

Next, we remark that 



^-V s \\ H - 1/2{2l) <Ce. (93) 



0x2 0x1 
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and on E, using Theorem 9, we obtain 

.aui 



Ik £ ||h-i/2(E) < 

A direct calculation shows that 



dx 2 



dxi 



[ dx~ 2 



h-i/2(e) +Ce < ||^^-|| H -i/2( S) + Cs. 
2 -^IHloc(e) <Ce, 



(94) 



and our result is valid for the traces taken from either unconfined side or from the side corresponding 
to the porous medium. □ 

4 Appendix 1: Very weak solutions to the Stokes system in 

Let Gi e L 2 (^i) 2 , G 2 € L 2 (£li) 4 , and £ e L 2 (S) 2 . We consider the following Stokes system in Oi: 

-Ab + VP = Gi + div G 2 in fi l5 
div b = in Oi, 
b = £ on E T = SU{i 2 = h}, 
, {b, P}, is L-periodic in xi. 



(95) 



Our aim is to show the existence of a very weak solution (b, P) <G L 2 (f2i) 2 x P" 1 (fii) to problem 
(95). To this end, we use the transposition method from [6]. 

Thus, let us test problem (95) with a smooth test function (3>,7r), satisfying $ = on St and 
being L-periodic in X\. Furthermore, ir is L-periodic in x\. We obtain 



< Gi + div G 2 , * >=< - div (Vb - PI), * >= - / P div * dx + 
/ (2D(&) - irl)v£ dSdt + f b ■ ( - A* + Vtt ) dx. 

Jn T Jill V / 



(96) 



Let (g,s) e TL r 9- 2 ^ r (r2i) 2 xVF' J - 1 ' r (r2i), 1 < r < +oo, 1 < g < 2, and U = {z e VT«- 1 ' r (fii), z 
0}, and denote by "H* its dual. Let now {<fr,7r} be given by 



' -A* + V7T = g in fli, 

div <& = s in Oi, 
, $ = 0, on St, {<l?,7r} is L-periodic in xi. 



(97) 



After the elliptic regularity for the Stokes system in [25], for q ^ 1+1/r we obtain <I> e VK 9 ' r (fii) 2 , ir £ 
W 9_1 ' r (ni), with J" n 7r = 0, and the following estimates hold 



l*llw>.>-(fil) 2 + ||V7r|| W <,-2,r(f 2) dt < C (||g|U<,-2,r (ni) 2 + || Vs|| W «-2,r (fi,)) . 



(98) 
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Now, analogously to the approach in [6] where the stationary Stokes system was treated, for 
q > 1 + 1/r, we consider the linear form 

*(g, s) = (Gi + div G 2 , *) ni - (V* - 7tI)u)j: t , (99) 

where (*,tt) is given by (97). Since (*,tt) satisfies (98), the linear form £ : W q - 2 > r {Q.{f x H -> R 
is continuous, and we set 

Definition 14. (04 very weak variational formulation for the Stokes problem (95)). {b, P} is a very 
weak solution of the problem (95) if 

{b, p} g w 2 -«> r /( r - 1 )(n 1 ) 2 x U* (100) 

and satisfies 

(g,b) ni -(P,8) n ., n =e(g,s), Vgei^VseH. (101) 

Because of the linearity and continuity of £, Riesz's theorem implies 

Proposition 15. Lef 1 < r < +oo and 1 + 1/r < g < 2 and < £2, 1 >s T = 0. Then, there exists a 
unique very weak solution {b,P} for (95). It satisfies the following estimates 

||b|| W 2-«,r/(r-l)( ni )2 < cj || Gl || L l ( Ql )2 + || G 2 || W l-9,r/(r- 1) )4 + ||^|| ^l + l/r-9,r/(r-l) ( Et )2 | . (102) 

Another approach is to use directly the result from the article [9], which reads 

Proposition 16. Let Gi = and G2 = 0. Then for £ e L 2 (St), J e C2 = 0, t/iere exis<s a unique 
very weak solution {b, P} 0/ satisfying the following estimates 

ll b llif 1 /2 ( n 1 )2 < c||^||z,2 (Et) 2. (103) 

Furthermore, 

|||a; 2 | 1 / 2 Vb|| i 2 (0l) 2 + |||a; 2 | 1 / 2 7r|| i 2 (0l) 2 < c||^|| L 2 (Et) 2. (104) 

5 Appendix 2: Navier's boundary layer and compressibility 
corrections 

In this Appendix, for completeness of the paper, we recall the derivation of Navier's boundary layer 
developed in [11] and [12] and presented also in [14]. 

As observed in hydrology, the phenomena relevant to the boundary occur in a thin layer sur- 
rounding the interface between a porous medium and a free flow. In this Appendix we are going to 
present a sketch of the construction of the main boundary layer, used for determining the coefficient 
a in (5) and the coefficient in the interface pressure jump law (24). Since the law by Beavers 
and Joseph is an example of the Navier slip condition, we call it Navier's boundary layer. 

In addition to the notations from subsection 2.1, we introduce the interface S = (0, 1) x {0}, the 
free fluid slab Z + = (0, 1) x (0, +oo) and the semi-infinite porous slab Z~ = U^ =1 (Y F - {0, k}). The 
flow region is then Z BL = Z+ U S U Z~ . 

We consider the following problem: 
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Find {/3 b ,<j bl } with square-integrable gradients satisfying 

-A y /3 bl + V v u u =0 inZ + \JZ- (105) 

div y p bl = in Z+ur (106) 

[l3 bl ] s (-,0) = and [{V y /3 bl - cu bl I}e 2 ] g (., 0) = e 1 on 5 (107) 

/3 W =0 on U^Li (dY" s - {0, k}), {/3 U , u hl } is 1 - periodic in yi (108) 

By Lax-Milgram's lemma, there exists a unique f3 bl £ L 2 oc (Zbl) 2 , V y z € L 2 (Zbl) 4 satisfying (105)- 
(108) and uj u e Lf oc (Z + U Z~), which is unique up to a constant and satisfying (105). We note that 
due to the incompressibility and the continuity of (3 bl on S, considering V(3 bl or the symmetrized 
gradient (V + V*)/? w is equivalent. 

The goal of this subsection is to show that system (105)-(108) describes a boundary layer, i.e. 
that /3 bl and u bl stabilize exponentially towards constants, when \y 2 \ — > oo. Since we are studying 
an incompressible flow, it is useful to prove properties of the conserved averages. 

Lemma 17. ([11]). Any solution {/3 bl ,uj bl } satisfies 



L 



f /3 bl ( Vl ,b) d yi =0, V6 e R and f uj b \y 1 ,b 1 ) d Vl = f oj bl ( yi ,b 2 ) d Vl , Wh >b 2 >0, 
o Jo Jo 



(109) 

f ffivuh) dm = ( p\ l {yiM) dyi = - [ \V(3 bl (y)\ 2 dy, V6i > b 2 > 0. (110) 

Jo Jo JZbl 

Proposition 18. ([11]). Let 

Cf = / ff(vi,0)d yi . (Ill) 
Jo 

Then, for every y 2 > and y x G (0,1), \/3 U (y 1 ,y 2 ) - (Cf,0)| < Ce~ Sy2 , V<5 < 2tt. (112) 
Corollary 19. ([11]). Let 

C b J= f LO bl (y u 0)d yi . (113) 
Jo 

Then, for every y 2 > and Vl e (0, 1), we /lave | uj bl ( yi ,y 2 ) - C bl |< e" 27 ^ 2 . (114) 

In the last step we study the decay of [3 bl and uj m in the semi-infinite porous slab Z~ . 

Proposition 20. (see [11], pages 411-412). Let (3 bl and uo bl be defined by (105)-(108). Then, there 
exist positive constants C and j n , such that 

\P bl (yi,y2)\ + \V/3 bl ( yi ,y2)\<Ce-^, for every ( yi ,y 2 ) € Z~ . (115) 

Furthermore, the limit = lim — \ — - [ uj b \y)dy exists and it holds 

fc^-oo I Y F | J Zk 

|w W (l/i,lte)-«oo| <Ce~^, for every (y u y 2 ) e Z~ . (116) 

Remark 21. Without loosing generality, we take — 0. If the geometry of Z~ is axially symmetric 
with respect to reflections around the axis y± = 1/2, then C^j = 0. For the proof we refer to 
[13]. In [13] a detailed numerical analysis of the problem (105)-(108) is given. Through numerical 
experiments it is shown that for a general geometry of Z~ , C% ^ 0. 
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It is important to be sure that the law by Beavers and Joseph does not depend on the position 
of the interface. We have the following result 



Lemma 22. Let a < and let (3 a ' bl be the solution of (105)-(108) with S replaced by S a — (0, l)x{a} ; 
Z+ by Z+ = (0, 1) x (a, +oo) and Z~ = Z BL \ (S a U Z+). Then, it holds 

c a,bl = C bl _ a (n7) 

This simple result implies the invariance of the obtained law on the position of the interface. It is in 
agreement with the law of Saffman for the slip coefficient formulated in [21]. The law was confirmed 
numerically by Sahraoui and Kaviany in [22]. For more discussion, we refer to the book [15], page 
74, formulas (2.193) - (2.195) and page 81, Fig. 2.22 and formula (2.211). 

The reminder of the section is devoted to auxiliary functions correcting the compressibility effects. 
We define Q w , by 

div y Q bl (y) = p bl (y) - C bl H(y 2 ) in Z+ U Z" (118) 
Q w = on UgLi (dY s - {0, fc}), Q bl is 1-periodic in y x (119) 

[Q bl ]s = e 2 [ (C bl H(y 2 )~f3 bl (y))dy = -e 2 [ /3 bl (y) dy. (120) 

J Zb l J Z 

Proposition 23. (see [11], page 411) Problem (118)-(120) has at least one solution Q bl e ^(Z+U 
Z-) 2 C\C% C {Z+\JZ-) 2 . Furthermore, Q w G W^{Z+) 2 , Q w e W^(Z-) 2 , for all q e [l,+oo) and 
there exists 70 > such that 

e 7os/ 3 QM € H 1 {Z + UZ~) 2 . (121) 
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